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SPECTRAL GEOMETRY OF FLAT TORI 
WITH RANDOM IMPURITIES 

HENRIK UEBERSCHAR 


Abstract. We discuss new results on the geometry of eigenfunctions in disor¬ 
dered systems. More precisely, we study tori R^/LZ^, d = 2, 3, with uniformly 
distributed Dirac masses. Whereas at the bottom of the spectrum eigenfunc¬ 
tions are known to be localized, we show that for sufficiently large eigenvalue 
there exist uniformly distributed eigenfunctions with positive probability. We 
also study the limit L ^ oo with a positive density of random Dirac masses, 
and deduce a lower polynomial bound for the localization length in terms of 
the eigenvalue for Poisson distributed Dirac masses on R^. Finally, we discuss 
some results on the breakdown of localization in random displacement models 
above a certain energy threshold. 


1. Introduction 

In 1900 Paul Drude introduced a classical model mis] for the motion of electrons 
in a material with the aim of studying transport properties such as conductivity 
of metals. In particular Ohm’s law could be derived from the model. In 1933 
the model was supplemented with results from quantum theory by Hans Bethe and 
Arnold Sommerfeld as what is now known as the Drude-Sommerfeld or free electron 
model. 

Anderson discovered in 1958 that at low energy electronic transport could break 
down in disordered media, provided the disorder is sufficiently strong [T]. In par¬ 
ticular the eigenstates in this regime are exponentially localized. This phenomenon 
is today known as Anderson localization. A key question concerns the transition 
from a localized to a delocalized regime when the disorder becomes small compared 
with the energy. In 1979, Abrahams, Anderson, Licciardello and Ramakrishnan 
proposed their scaling theory which suggests that such a transition should exist in 
dimension d > 3, the case d = 2 being critical, whereas for d = 1 there is only 
a localized regime. Regarding d = 2 the widely held belief is that no transition 
exists, however the localization length can be very large compared with the size of 
the system. 

The type of disordered system considered in this paper is a box A = [L/2, L/2]'^, 
L > 0, which contains N independently uniformly distributed impurities, which are 
modeled by Dirac masses (also known as Fermi’s pseudo-potential, delta potentials 
or point scatterers). In this paper we will consider periodic boundary conditions, 
i.e. we will study fiat tori However, our results can easily be 

generalized for Dirichlet or Neumann boundary conditions. 
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2. Tori with random impurities 

By a random Schrodinger operator on we mean the following type of stochas¬ 
tic differential operator 

( 2 . 1 ) Hn = -^+Y,V{x-uj), V&C\Ti), 

CJ G 

where fl, |0| = TV, is a stochastic process on T^. 

A simplihed model H replaces V = 5, where i5 denotes a Dirac mass cen¬ 
tered at the origin. Such a “potential” is known as Fermi’s pseudo-potential, a 
delta potential or a point scatterer. The model approximates the generic random 
Schrodinger operator iJn well in the regime, where diam supp V <C w{E ), and where 
w{E) X Ij^fE denotes the wavelength, and E the energy. 

The advantage of considering the formal operator H is that it can be realized 
as a self-adjoint extension of the restricted Laplacian Hq = — A|(^oo(-.]rd _q), which is 
a positive symmetric operator with deficiency indices {N,N). The extensions are 
parametrized by a unitary matrix U G U{N), and the corresponding self-adjoint 
extension is given by the restriction of the adjoint E[q to the domain of functions 
/ € E[^ (T^ — D) which are of the form 

fix) = 9{x) + (v,Gi(x)) -b ([/T,G_i(a;)) 

and G±i(x) = (G±i(a;, Wi), • • • , G±iix, lom)), where Gi{x,uj) and G_i(a;,a;) are in¬ 
coming and outgoing circular waves at the point w € T^. 

We note that the parameter space of self-adjoint extensions is much larger than 
the physical parameter space which is of dimension N (there are N real coupling 
constants). In fact a non-diagonal matrix U corresponds to an extension which 
violates local conservation of mass in the scattering process (i. e. it is possible for 
part of the wave to enter into one impurity and emerge from another, as if a wire 
were attached between the two points). Assuming local conservation of mass in 
all scattering processes we reduce the parameter space to the subgroup of diagonal 
unitary matrices D{N) C U{N). Note that the operator E[ corresponds to a matrix 
D = e“^Id, for some (p € (—7r,7r), since all coupling constants are equal to 1. 

3. Anderson localization 

Let us consider a particular type of random Schrodinger operator on M’’*, a ran¬ 
dom displacement model (cf. Figure 1) 

(3.1) _A-f ^ D(x-C-W 5), I/eGr(M‘’*) 

eez-i 

where the displacements wj are i.i.d. random variables with radially symmetric 
probability density P(x) = T(|x|), P(0) > 0 and P € G“(M+) strictly decreasing. 
The disorder in the present system can for instance be measured by the average 
displacement, Sq = Ep(Al), where A is a random variable with probability density 

P. 

If the probability distribution P were given by a Dirac distribution centered at 
the origin, i.e. if Z'’* -b a; = and So = 0, then we know from Bloch’s theorem 
that the spectrum has a continuous band structure with corresponding extended 
eigenstates. However, if there is sufficiently strong disorder, i.e. (5o is large enough. 
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Figure 1 . The picture above shows a random displacement model 
in two dimensions. In black the lattice in grey the randomly 
perturbed lattice I? +uj. 


then we have, almost surely, pure point spectrum and the corresponding eigenstates 
are exponentially localized at the bottom of the spectrum. 

Definition 3.1. Let E 2 > Ei > 0. We say that Hq is exponentially localized on 
[Ei,E 2 ], if Hq has a.s. pure point spectrum on [Ei,E 2 ] and the eigenfunctions ifx 
satisfy the bound : 

Bta G : Vx G 

where the localization length Lioc depends on Ei , E 2 and the localization centre xx 
as well as the constant C depend on the sample w G if. 

Mathematical proofs of Anderson localization were for instance obtained by 
Goldsheid, Molchanov and Pastur Frohlich and Spencer [5], Simon Bour- 
gain and Kenig [3], Klopp m and Germinet [7]. However, a crucial question 
concerns the existence of a delocalization transition in the regime where the disor¬ 
der is small compared with the energy, in the case of a random displacement model: 
(5o < E. 

In 1979, the “Gang of Four”, Abrahams, Anderson, Licciardello and Ramakrish- 
nan [2] proposed what is now known as the “scaling theory of localization” which 
predicted that the answer should depend on the dimension d. For d = 1 the the¬ 
ory predicts that there is only a localized regime, i.e. any amount of disorder is 
sufficient to localize all eigenstates. For d > 3 it predicts the existence of a delo¬ 
calization transition, i.e. when the disorder is sufficiently small compared with the 
energy, the eigenstates should become delocalized. The case d = 2 was identified 
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as critical, and it was famously conjectured that there should only be a localized 
regime, however, the localization length could be very large. 

The mathematical characteristics of the delocalized regime are continuous band 
structure above a certain energy threshold and associated extended eigenstates. A 
very interesting phenomenon is observed near the transition point between the two 
regimes, where physicists have observed a multifractal structure of the eigenstates 
0 . 

Very few rigorous mathematical results exist regarding the multifractal or de¬ 
localized regime. We will continue to describe some rigorous results about the 
geometry of the eigenstates of Poisson random Schrodinger operators with Dirac 
masses on flat tori. In particular these compact models allow us to give polynomial 
lower bounds on the localization length on In fact for certain stochastic pro¬ 
cesses, such as random displacement models, one can go much further and prove 
the breakdown of exponential localization for sufficiently large eigenvalues [T3] . We 
will sketch these results at the end of this paper. 

4. Spectrum and eigenfunctions 

We consider a flat torus T^, d = 2,3 with N independent uniformly distributed 
delta potentials located at the points D = {xi, ■ ■ ■ , XArjand we let U = e“^ Wat. We 
denote the associated self-adjoint extension of by —A,^. We denote 

the Green’s function, the resolvent kernel of the Laplacian on by 

Gx{x,y) = ^ S{x - y), A ^ cr(-A). 

The spectrum of —A,^ consists of two parts, old eigenfunctions of the Laplacian 
which vanish on D and new eigenfunctions which, almost surely, diverge at each 
point in D. The divergence is of order log \x — Xj \ in d = 2 and of order l/|a; — Xj\ 
in d = 3. 

The new eigenvalues are solutions of the equation 

(4.1) detA^ = 0 
where the matrix entries are given by 

= Gx{xk,xi) - ^Gi{xk,xi) - tan(^)3Gi(a;/c, xi). 

In particular one can show that per old Laplacian eigenspace the self-adjoint exten¬ 
sion —A^ produces at most N new eigenvalues, which almost surely lie in between 
two neighbouring Laplacian eigenvalues. 

The associated new eigenfunctions are superpositions of Green’s functions 

(4.2) tpxix) = {v,Gxix)), iiGkerA^ 

and we recall (GA(a;) = {Gx{x,xi),- ■ ■ ,Gxix,XN))- Also note that u is a function 
of the random variables Xi, • • • ,xn. Given det = 0, the matrix has almost 
surely rank N — 1 which implies dimker A^ = 1, so the choice of v is unique up to 
normalization. 


5. Uniformly distributed eigenfunctions 

The aim of this section is to given information on the geometry of the new 
eigenfunctions of the random operator — A,^. For a large number N of impurities 
we expect the eigenfunctions to be exponentially localized at the bottom of the 
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spectrum. However, we will show that for sufficiently large eigenvalues there exist, 
with a certain positive probability (which tends to zero as N ^ oo), eigenfunctions 
which are uniformly distributed on the torus T^. In particular this fact gives 
information on the size of the localization length, yielding a certain polynomial 
lower bound. 

Before we state our main result we recall that each new eigenvalue A satisfies 
almost surely A € (Aj, Aj+i) for some j. We have the following result. 

Theorem 5.1. Denote the Laplacian eigenvalues on by Fix any gq € 

There exists Jd = {jfc}fe C N, an index subset of full density such that 
for any j G Jd and any new eigenvalue of —A G (A,-,A,+i), we have with 




We may now apply this result to give a lower bound on the localization length. 
The equidistribution theorem above implies that if we are still in the localized 
regime, then the localization length must exceed the size of the torus: 


-b/oc -b. 


Let us now fix a positive density of impurities 

N = pL^ 


and observe that our stochastic process on converges to a Poisson process of 


density p on in the limit L —>■ oo. Our equidistribution result on the torus of 


therefore gives information on the localization length for a Schrodinger operator 
with Poisson delta potentials on 

In particular, equidistribution occurs when 




where 



And this implies that the localization length must satisfy the lower bound 

Lioc> A“L 


6. Elements of the proof 


Consider the square torus = K^/27rZ^. 
Let 


N 




We introduce a test function a G and the L^-normalized eienfunction is 

denoted by tpx = 4'a/||^'a||2- 


^By which we mean 'iC, ■. |a(C)l < l“o(C)l- 
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We have 





were 

ec(a^) = 

Our goal is to show that a matrix element correponding to a nonzero mode 
C 7 ^ 0i is small if A > 1. 

Let us define the annulus 


Ac(A, L) = U G Z 2 I lie - Cp - A| < L}, C e ^ 2 , L > 0 . 


Using Chebyshev’s inequality one can show that with probability > ^ we have 


I (e^-0A,'*/'A) P < ^ 


S?eAlo(A,A'*)(l^l^ “ 



where i5 G | — 0) and 0 = |y| is the best known exponent (due to Huxley [5]) 
in the circle law 


N{X) = #{|e|^ < X I e G Z^} = ttX + 0,(X®+"). 


For a density one subsequence of new eigenvalues A we have that ^ € Ho(A, A"^) 
implies 

e^Hc(A,A^)o||^-Cp-A| > A^ 

so that 


ll?-CP-A| 

Also the error term in the circle law implies the following bound on the number of 
lattice points in a thin annulus 


#Ao(A,L)=0,(A«+^). 


Combining the bounds above we obtain the estimate (recall S > ^9) 

(ie-CP-A)-2 = 0.(A-2'5+«+^). 

||5|2-A|<A^ 

Furthermore, we have for a density one subsequence of eigenvalues A the lower 
bound 


E (ICp-A)-2>C(e)A-C 

5GAo(A,A'5) 

So, for generic A, the following bound holds with probability > 

l(ecV'A,i/'A)PS iVA-i+3«+C 


j_ 

N 


This argument can easily be extended for any trigonometric polynomial with 
nonzero frequencies. The result for a general torus T| = K^/LZ^ then follows by a 
simple scaling argument. 
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7. Delocalization for random displacement models 


Let B[^ = [—L,L]‘^ and consider 
Hq.l = —A 


{GZ'^nSj 




X > 0 


with Dirichlet boundary conditions, where the displacements uj^ are i.i.d. random 
variables with probability density P G C°(K‘^) and suppP C B(0, |). 

Denote by ipx an L^-normalized eigenfunction of Hq ^. Fix y G (7° 
and llxlli = 1, and introduce the smoothed eigenfunction 

1/2 

X{x' - x)\^j;{x')\^dx' 

'Bl 




One of the problems with Definition 3.1 is that the localization centre depends on 
the sample w G fl. To observe localization it is therefore advisable to consider the 
two-point correlation function of the smoothed eigenfunction which is dehned 
for any two points x,y G B^ by 

Qx{x,y) = '^/{{x)'S^{y). 

The decay is then observed with respect to the distance \x — y\, which is independent 
of the sample w. 

We have the following alternative definition of localization. 


Definition 7.1. Let F G C'°(K+) be strictly decreasing and 
||'0aII2 = 1- say that is F-localized on an interval I C K+ if, for suf¬ 

ficiently large L, we have 

'dx,yGBL-. E < F{\x - y\) 

and the limit of the LHS as L —>■ oo exists. 

We can show that for large enough energy fails to be localized. We denote 
the set of distinct eigenvalues of the Dirichlet Laplacian on 73^ by A^. 

Theorem 7.2 (H.U. 2015). There exists Eq ^ 1 and a full density subsequence 
Af C Al such that for any Xk G Af, Xk > Eq, 


El ^ 0^(21,2/) >1. 

A'^C(Afc,Afc+i) J 

In particular, this result implies that the localization length blows up at a certain 
critical energy threshold: Ljoc oo as E ^ Eq. 

7.1. The limit of large tori. The key idea in proving the above theorem is to 
study the operator iJn.i in the limit L —>■ oo. One of the key obstructions to doing 
this is the dependence on the number of potentials, N xi L‘^, in the error term in 


equation (5.11. 


However, we are able to improve the estimate of the error term significantly in 
the case of random displacement models: Fix oq G C°°{Bi) and any e > 0. Denote 
by A = {Ajl^g the set of distinct eigenvalues of the Dirichlet Laplacian on Bi. 
There exists a full density subsequence A' C A such that we have, for sufficiently 
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large Afe G A'^, with probability 1 — e for any a € |o| < |ao| and any 

A G (Afe, Afe_|_i) 

^ b{y)\^{{y)\^dy=^^i^j^ b{y)dy + , b{-) = L-'^a{-/L). 

Now let Oq = X G ||x||i = 1 and suppx C i?(0, eg) for some small eg. 

We thus obtain the lower bound 

for any L 3> 1. In particular, we have the following lower bound for the two-point 

correlation function 

<dx{x,y)>\x-y\~'^, \x-y\xiL. 

In particular if we sum over all A G X = (Afc,Afc+i), Afc G A^, and take the 

expectation, we may show 

E(^e^(x,y)) >1 

\Aei / 

where we have used that 

#{AgX}x 
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